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Abstract. We consider in C n with coordinates £ = (£l ,-..,£«) the class of symmetric homogeneous qua- 
dratic dynamical systems. We introduce the notion of algebraic integrability for this class. 

We present a class of symmetric quadratic dynamical systems that are algebraically integrable by the set 

of functions hi[t), , h n (t) where hi(t) is any solution of an ordinary differential equation of order n and 

/ifc(t) are differential polynomials in hi(t), k = 2, . . . , n. We describe a method of constructing this ordinary 
differential equation. We give a classification of symmetric quadratic dynamical systems and describe the 
maximal subgroup in GL(n, C) that acts on this systems. 

We apply our results to analysis of classical systems of Lotka-Volterra type and Darboux-Halphen system 
and their modern generalizations. 

Introduction. 

Our study of polynomial dynamical systems (see ]T|) drew our attention to the class of symmetric quadratic 
dynamical systems (see definition below) . For such systems we introduce the notion of algebraic integrability 
and in the generic and almost generic case find solutions of the algebraic integrability problem. Our approach 
reduces the problem of integrability of a symmetric quadratic dynamical system to the question of solving 
an ordinary differential equation. This opens the way to apply ideas of Kovalevskaya and Painleve: to 
find solutions of dynamical systems in terms of special functions which are solutions of remarkable ordinary 
differential equations. The classical example here is the solution of Darboux-Halphen system in terms of 
Chazy equation solutions. 

In the first part of the article we consider for any n the space of quadratic dynamical systems and the 
canonical action of GL(n, C) on it. We describe the subspace of symmetric quadratic dynamical systems and 
the maximal subgroup in GL(n, C) acting on the space of such systems. In the second part we introduce the 
notion of algebraic integrability and give an algorithm of solving the problem of algebraic integrability for 
generic and almost generic symmetric quadratic dynamical systems. This algorithm reduces the problem of 
algebraic integrability to the problem of solving an ordinary differential equation. We give the action on such 
equations of the normalizer subgroup in GL(n, C) of generic symmetric quadratic dynamical systems. In the 
third part we show how our method works for known problems of quadratic dynamical systems theory. 

The authors are grateful to E. B. Vinberg for useful discussions. 

1. Symmetric quadratic dynamical systems. 

1.1. Quadratic dynamical systems. 

Fix some n£N. Let us consider in C™ with coordinates £ = (£i, . . . , £„) T the general quadratic dynamical 
system 

(1) &(<)=4 J '&(W). A? = A?* = const, k = l,...,n. 

Here and below we use the usual Einstein summation convention. We identify the space of quadratic 
dynamical systems with the \n 2 (n + l)-dimensional linear space of tensors A = (A^ 3 ) such that A 1 ^ 3 = A 3 ^ 1 . 
Note that the system ([T]) is homogeneous with respect to the grading degi = 4, deg£i = —4, degA^ 3 = 0. 
The system ([T]) defines a linear operator C = C{A) : C[£i, ...,£«] — > C[£i, . . . , £„] with deg C — — 4 as 

d 
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Note that P(£) is an integral of system {jXJ) if and only if £P(£) = 0. 

The group GL(n, C) acts by linear change of coordinates i] = B£ where B — (Bj) E GL(n, C) on the space 
of dynamical systems (TJJ. This action is a representation of GL(n,C) on the space of tensors A = (AV): 

A? ^AYBl{B-%{B^y r . 

1.2. Symmetric quadratic dynamical systems. 

In this section we introduce three definitions of a symmetric quadratic dynamical system and prove their 
equivalence. 

The symmetric group S„ is the group of permutations of n elements. We identify S„ with the subgroup in 
GL(n,C) generated by the matrices obtained by a transposition of two lines in the identity matrix E, thus 
S n acts on C" by permutation of coordinates (£1, . . . , £ ra ). 

Definition 1. A system (flj is called symmetric if S n lies in the stabilizer of this system. 

Definition 2. A linear operator L : C[£i, . . . ,£ n ] -> C[£i, . . . ,£„] is called symmetric if T(LP(£)) = LTP(£) 
for any T E S„ and any polynomial P(£), where TP(£) = P(T (£_)). 

Definition 3. A system (JlJ is called symmetric if C(A) is symmetric. 

Let Sym C C[£i, . . . , £ n ] denote the ring of symmetric polynomials. According to classical results, Sym is 
isomorphic to a polynomial ring with n generators. The best known generators of this ring are the Newton 
polynomials and the elementary symmetric functions. We denote the Newton polynomials as pk — Pk{£) = 
and the elementary symmetric functions as <Jk = Cfc(£) = V ( CuC^ • • -£ih- ^ or simplicity set 

Po = n. One can pass from one of this sets of generators to another by the classical identities (see [2]) 
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The universal algebraic map 
(2) S:C"->C n : £.->(*!,... on) 

is the projection to the orbit space of the action of S n on C™. We have Sym 
Definition 4. A system |T]) is called symmetric if CP E Sym for any P E Sym. 
Theorem 5. The definitions [TJ [3] and [4] are equivalent. 

Proof. For T E S„ denote by T(z) the image of i under the permutation T. We have 

d_ 

and 

T£P(0 = (T£)(TP(0), 

thus definition Q] is equivalent to T£ = which is equivalent to defmition[3J Let system ([T]) satisfy definition 
[3] Any T e S„ does not change symmetric polynomials, moreover, this is a characteristic property of such 
polynomials. Each permutation T does not change the operator C(A). Hence for a symmetric polynomial 
P the polynomial CP is symmetric. 

Let system ([!]) satisfy definition |4j The polynomial CP = A^ J £,i£,j(§^P is symmetric for any symmetric 



(r£)PK)=4j| T(i) ^^p(e), 



polynomial P, and thus does not change under the action of S„. Therefore, we have A^£j F siCj'afjP = 



Ai j ^-±-P for any T E S„ and any symmetric P. We conclude that (a^™ ~ Ci",-^/' = for 
any symmetric P. Now if for a derivation L we have LP = for any symmetric P, then S'*L = 0. But S is 



a local diffeomorphism, hence A 



T(i),T(i) 
T(k) 



A! J for any T E S„, and we come to definition [TJ 



□ 



Lemma 6. Each symmetric quadratic dynamical system has the form 
(3) Ck(t) = ^ j ^(tMt), fc=l,...,n, 

where A™ = A" for any i,j,k, A™' m = A™!'" 1 ', A™' TO = A™?' m ' , A%>» = A™/™' and A*'" = A™ V /or 
any n, m, I, n' , m' , I' such that n ^ m, I ^ m, I ^ n, n' =fi m' , I' ^ ml , V 7^ n! . 

Proof. Let T 6 S n be a transposition of i and j. In the general quadratic dynamical system ([T]) according to 
definition [D we have Af = Af , Ajf = Atf, A\' 1 = Af , A\> k = Af \ A k ' 1 = A^ for any k ^ i, k ? j, I ± i, 
I ^ j. By applying this equalities for all i, j we see that the only possibly pairwise different coefficients of 
the system are A™' ro , A%> m A™'™ and A™' 71 for n ^ m, I ^ m, n. □ 

Corollary 7. Each symmetric quadratic dynamical system has the form 



. n. 



(4) a = ^i+KkPi(o+iPi(o 2 +sp2(o, k=i, 

for some coefficients a, f3, 7, 5. 

Proof. The form (0| is a linear change of variables of ([3]), namely A™ : ™ 1 = a + (3 + 7 + S, A" l,m = 7 + S, 
2A™> n = (3 + 2 7 , A\ n ' n = 7. ^ □ 

Lemma 8. The stabilizer subgroup in GL(n, C) of symmetric quadratic dynamical systems coincides with 
S n • 

Proof. Consider the system 

€k = Cfei k = 1, . . . ,n. 

Each g in the stabilizer subgroup does not change this system. Let g bring £/. into rjk = a^Cs- The dynamical 
system in r]k is 

Vk=Vk, k = l,...,n. 

Thus for any £ we have 

i<j 

therefore a\ ^ only for one value of i for each k and a\ = (a l k ) 2 . Using g G GL(n, C) we conclude that 
there exists a permutation T 6 S„ such that a,k,i = 1 if i = T(k) and dk,i = otherwise. □ 

Definition 9. The system ([I]) is said to be quasi- symmetric with respect to B £ GL(n, C) if this system in 
the coordinates n = P£ is symmetric. 

1.3. Integrals of symmetric quadratic dynamical systems. 

Lemma 10. Let P(£) be an integral of a symmetric quadratic dynamical system (J3]). Then the polynomials 
T(P(£)) = P(T(£)) where T € S n are integrals of the same system. 

Proof. According to definition [5] we have £(TP(£)) = T(£P(£)) = 0, thus TP(£) is an integral of the 
symmetric system (J3J). □ 

Remark 11. Pac/i integral P(£) 0/ a symmetric quadratic dynamical system ([3]) can fee presented as the sum 
of its graded components, one component in each grading, each component being an integral of the system. 



Let us consider the projector ir : C[£i, . . . , £ n ] — > Sym defined as 

= h E Tp («- 



n 
tss„ 

Corollary 12. Lef P(£) oe a homogeneous integral of a symmetric quadratic dynamical system ([3]). XTien 
£/ie polynomial P(£) = ttP(£) * s a homogeneous symmetric integral of the same system. 

Proof. Obviously, a sum of integrals of a dynamical system is itself an integral of the dynamical system. 
Therefore the corollary follows from Lemma ITOl □ 



1.4. Generic symmetric quadratic dynamical systems. 

In this section we introduce the notion of generic symmetric quadratic dynamical systems and classify 
such systems for every n. 

For homogeneous generators a±, . . . , a n of the ring Sym with deg au — — 4fc a symmetric system @ implies 
a homogeneous dynamical system 

(5) a' k — Ca k , k = l,...,n. 
Using the grading one can rewrite ([5]) in the form 

(6) a' k = 5fe+i(oi,...,Ofc) +Cfc<Zfc+i, fc = 1, ...,n, 

where c n = and gj{a%, . . . , Oj-i) are homogeneous polynomials with deg <7j = — 4j. 

Definition 13. A homogeneous symmetric system ^ is generic if Cfc ^ for fc = 1, . . . , n — 1. 

Lemma 14. j4 system being generic does not depend on the choice of generators in the ring. 

Proof. Let a\,...,a n and b\, . . . , b n be two generators in the ring with deg dj = deg bj = — 4j. Consider the 
ideal J generated by decomposable elements of Sym. We have ak — Ckbk mod J, where Ck ^ 0, Ck S C. 
Thus, ([5]) implies 

o' fe = Gfe + i(6i, . . . , 6 fc ) + Ck—^-bk+i, k=l,...,n, 
for some homogeneous polynomials Gj(bi, . . . , bk). Therefore Ck i= c fc -77^ ^ 0. □ 
Theorem 15. For n = 1 each generic symmetric quadratic dynamical system has the form 
(7) 

For n = 2 eac/i generic symmetric quadratic dynamical system has the form 

(8) £fc=<fe+/?6pi(0+7Pi(0 2 , fc = l,2, 
wii/i a =/= 0. 

For n ^ 3 eac/i generic symmetric quadratic dynamical system has the form 

(9) £fc = «4 2 + ^Pi(0+7Pi(C) 2 +*P2(0. fc = l,...,n, 
uraift a^O, a + n<5 7^ 0. 

Proof. For n = 1 any quadratic dynamical system has the form <[7j) , it is symmetric and the condition for it 
being generic gives no restrictions. 

For n = 2 using corollary [7J and the relation P2 — Pi + 2£,iPi — 2^ = we get the form © . Therefore 

p 1 (C)' = (/3 + 2 7 )p 1 (e) 2 +ap 2 (e), 

which is the first equation of (JB]) for a./. — Pk- Thus is generic for a ^ 0. 
For n)3we get from corollary [7] directly. Therefore 

^ Pk (0' = otpk+iiO + PPk(OPi(0 + 7Pfe-i(0Pi(C) 2 + *U*-i(0P2(0, fc = 1,2, . . ., 
where Po(£) = n - Thus in ([H]) for = pk we have Ci = a + nS and Ck = a for fc = 2, . . . , n — 1. □ 



Definition 16. A quasi-symmetric with respect to B system is generic if the symmetric system in the 
coordinates rj = Bt; is generic. 
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1.5. Almost generic symmetric quadratic dynamical systems. 

Definition 17. A homogeneous symmetric system (|3]) is almost generic if in ([6]) we have Ck ^ for 
fc = l,...,n-2. 

Lemma 18. A system being almost generic does not depend on the choice of generators in the ring. 

The proof coincides with the proof for Lemma [HI 

Theorem 19. For n — 2 each symmetric quadratic dynamical system is almost generic. 

For n = 3 each almost generic symmetric quadratic dynamical system has the form ^ with a + nS ^ 0. 
For n 4 each almost generic symmetric quadratic dynamical system is generic. 

The proof is the proof of Theorem 1 1 5 1 with obvious modifications. 

Definition 20. A quasi-symmetric with respect to B system is almost generic if the symmetric system in 
the coordinates r\ = B£ is almost generic. 

1.6. Action of GL(n,C) on symmetric quadratic dynamical systems. 

Recall E is the n x n identity matrix and denote by e the n-dimensional vector e = (l,...,l) T . 
Consider the set of matrices B(X, q) = XE + qee T . With respect to the standard multiplication of matrices 
it is a commutative monoid with multiplication 

B(Ai, gi)B(A 2 , q 2 ) = B(AiA 2 , X ± q 2 + A 2 ?i + nqiq 2 ). 

We have det B(X,q) = A™" 1 (A + nq), therefore B(X,q) £ GL(ra,C) if A ^ 0, A ^ -nq. Denote by B the 
subgroup of matrices B(X, q) £ GL(n, C). The inverse element is 

Lemma 21. The centralizer ofS n in GL(n, C) is B. 

Proof. From the definition of B we see that B commutes with S n - Let G £ GL(n, C) satisfy T~ 1 GT = G for 
any T £ S n . Then G\ = G], G\ = G) for any i, j by the action of T = (i, j), and G\ = G k , G) = G{ for 
any i ^ j, i ^ k by the action of T = (j, k). Therefore G = XE + qee T for A + q = G\ and q — G\ for any 
i^j. " □ 

The following lemma describes the action of B on symmetric quadratic dynamical systems. 

Lemma 22. The action r] = £>(A,q)£ with A ^ 0, A + nq ^ brings each generic symmetric quadratic 
dynamical system (0| into the generic symmetric quadratic dynamical system 



a 9 2qa — X/3 q 2 (a + nd) — XfX'y — 2qS) 7 q(a + n8) + X5 
A % - X(xTn-qf kPM XHxTn^) + A^ 1 



Proof. We have n = A£ + qpi(t;) e i thus Piiv) = (A + ng)pi(C) ancl P2(»7) = A 2 p 2 (£) + (2A + nq)qp 1 (£) 2 . 
Substituting this relations into (j4|) we get the required system, which is generic for a^O and a + nS =/= 0. □ 

Corollary 23. The action r\ — B(\,q)t; with A ^ 0, A + nq ^ brings each generic symmetric quadratic 
dynamical system with n = 1 into the generic symmetric quadratic dynamical system 

l a 2 

and with n = 2 into the generic symmetric quadratic dynamical system 

{X + 2g)a M(\ + q )a-X 2 f3 g(X + q)a + X 2 1 2 

Vk - X 2 Vk X 2 (X + 2q) 71kPl ^> + X 2 (X + 2q) MV) ' fc_M ' 

Lemma 24. The normalizer subgroup in GL(n, C) of generic symmetric quadratic dynamical systems is the 
normalizer of S n in GL(n, C). 

Proof. According to lemma |5] the stabilizer subgroup is S„. Let h £ S n and let g lie in the normalizer 
subgroup in GL(n,C) of generic symmetric quadratic dynamical systems. Then h' = ghg~ x lies in the 
stabilizer subgroup S„. □ 
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Theorem 25. The normalizer subgroup in GL(n, C) of generic symmetric quadratic dynamical systems is 
equal to B x S„ . 

Proof. Let us denote the normalizer subgroup in GL(n, C) of generic symmetric quadratic dynamical systems 
byG. 

Form Lemma 1241 it follows that any element g G G induces an automorphism tp g (h) — ghg^ 1 of the 
group S„. It is well known that for n / 6 the group S n has no outer automorphisms, thus ip g is an inner 
automorphism. For the proof that ip g is an inner automorphism for n = 6 see below. 

Take any g G G. Since ip g is an inner automorphism, there is an element hi G S n such that tp g = 
(fh', i.e. ghg^ 1 = h'h(h')^ 1 for any h G S„. This can be rewritten as ((h')~ 1 g)h — h((h')~ 1 g). Therefore 
(h')~ 1 g G GL(n,C) commutes with any h G S n . By Lemma [2~T1 we get (h')~ 1 g G B. 

Thus G is generated by B and S n in GL(n, C). From B D S„ = {e} and bh = hb for any b G B and h G S„ 
we get G= (B,Sn) = B x S n . 

For n — 6 let ^ be a non-inner automorphism. This automorphism takes a transposition T± = (i,j) 
to a product of transpositions T 2 — (ii, i2)(*3j ii)(ih, ie)- It follows from the fact that ip preserves order 
and sign, but does not preserve cyclic type of permutation as a non-inner automorphism. Consider T\ 
and T2 as elements of GL(n, C). The element T\ has eigenvalues (i, —i, 1, 1, 1, 1), the element T 2 has eigen- 
values (i, —i, i, —i, i, —i). This implies that T\ and T 2 are not conjugate in GL(n, C). Hence, tp g is an inner 
automorphism of S„ as well for n = 6 and any g G G. □ 

Lemma 26. The orbits of the action of £>(A,q) on generic symmetric quadratic dynamical systems are as 
follows: 

For n = 1 there are two orbits: a ^ and a = in (|7|). 

For n — 2 and a ^ —/3 there is a one-parametric space of orbits with exactly one representative of the 
form 

(10) ti=e k +7Pi(0 2 , k = l,2, 

in each orbit, for a 7^ 47 there is a one-parametric space of orbits with exactly one representative of the form 

(11) & = $+MkPi(Z), k = l,2, 
in each orbit, for a — —j3 = 47 a representative is 

(12) z' k =&-tkPm+\pi(£) 2 , k=i,2. 

For n ^ 3 there is a two-parametric space of orbits with exactly one representative of the form 

(13) H' k = e k + kkPi(0+JPi(0 2 , k = l,...,n, 
in each orbit. 

Proof. For n = 1 this result follows from corollary |2"51 

For n = 2 for a ^ —f3 the non-degenerate matrix B(\,q) with A = \/ a(a + f3), q = \ A ^ A ~°^ , brings 
system ((8]) into the generic symmetric quadratic dynamical system (fTQ| with 7 = j ^a+p • ^ or a ^ ^1 the non- 



degenerate matrix B(X, q) with A = y/ a (a — 47), q = g — a , brings system ((5]) into the generic symmetric 
quadratic dynamical system (TTTI) with /3 = • For a = — [3 = 47 the non-degenerate matrices £>(A,q) 
with A ^ 0, q = I a ^ A q , bring system (|5J| into the generic symmetric quadratic dynamical system (fT 

a 2 (a-t-n5) 



For n ^ 3 the non-degenerate matrix B(\,q) with A = S 5 q = — -27° , Sa l , a ^ 0, a + ni^ brings 



system ([9]) into the generic symmetric quadratic dynamical system 

_. 2 a(2u5a - 26a 2 + (3aa + f3anS) , (5 2 a 2 - 5 2 a 2 + ja 2 a + •ya 2 n5) _, 2 ~ 

4 = "4 + — — ~. £kPi(0 + — — ~ Pt(Q + Sps(0> 

a z (a + nd) a^(a + no) 

k = 1, . . . , n. Thus in each orbit one can find a system with 5 = 1 and 5 = 0. We get the system ([j"3"|) 
with (3= ^2±M,7 = (* 2 +w»*). □ 



2. Algebraic integrability. 



2.1. Algebraic integrability of symmetric quadratic dynamical systems. 

In this section we introduce the notion of algebraic integrability. 
Consider the equation 

(w) ?-h 1 e- i + ... + (-i) n h n = o 

with the roots (£i, . . . , £„). By Vieta's formulas we have hk — Ofc(£)> k — 1, . . . , n, therefore (hi, . . . , h n ) — 
■ ■ ■ ,£n), where S is the universal algebraic map ©. The discriminant of (fT4")) is A„ = finite — fj) 2 - 
It is a symmetric polynomial of , £n, we have deg A„ = — 4n(n — 1). Denote by S the restriction of 

the universal algebraic map on {£: A„ ^ 0}. The map S is a covering. Therefore for each set of functions 
(hi(t), . . . , h n (t)) C ImS the set of functions (£i(t), • • • , £n(t)) C S'~ 1 (/ii(t), . . . , h n (t)) is uniquely defined 
by initial data. 

Lemma 27. for a symmetric quadratic dynamical system of the form @ one has 

£A„ = (n - 1) (2a + n/3) cr^A^ 

Proof. We have 

£A» _ v 2£(& - 6) _ V V 2(og + gfagi(fl + 7Pi(0 2 + M gfe - _ 

= 2 2^ , / =2 2^(ate + Ci)+/3pi(0) = (n - 1) (2a + nf3)p 1 ($). 

The remark Pi(£) = o"i(£) finishes the proof. □ 

This lemma implies each symmetric quadratic dynamical system defines a vector field C tangent to the 
discriminant curve {£ : A„ = 0}. Therefore for each symmetric quadratic dynamical system its restriction 
on the space {£ : A„ ^ 0} is well-defined. 

Definition 28. We say that a symmetric system ([3]) with initial initial data (£i(to)i • • • j£n(to)) 4- ^« is 
algebraically integrable by a set of functions (h%(t), . . . , h n (t)) if S(£\(to), . . . , £n(*o)) = (hi(to), . . . , h n (to)), 
the set (hi(t), . . . , h n (t)) lies in the image of S for any t and the uniquely- defined set (£i(t), • • • , £n(t)) C 
S^ 1 (hi(t), . . . , h n (t)) is a solution of the dynamical system ([3]). 

Problem 29 (Algebraic integrability) . Given a symmetric system of the form ([3]) , find an ordinary differen- 
tial equation of order n and differential polynomials hi, . . . , h n in hi such that for the each solution h\ of the 
ordinary differential equation each set of functions h±, h 2 , ■ ■ ■ , h n such that A(/i l7 . . . , h n ) ^ algebraically 
integrates this system. 

Theorem 30. For each generic symmetric system ((3|) with the initial data £(to) — (£i(io)) • • • >£n(to)) such 
that £i(to) ij(to) for any i ^ j there is a solution of the problem of algebraic integrability. 

Proof. Consider system ([6]) with = being the basis of elementary symmetric functions. Under the 
conditions of the Theorem we have Ck 0, k = 1, ...,n— 1. Hence, o~j(t) with j — 2,...,n can be 
expressed as polynomials in o~i(t), . . . , o~j-i(t) and their derivatives from the (j — l)-th equation. Thus, the 
last equation gives a single homogeneous differential equation for a\(t) in which the coefficient multiplying 
the highest derivative of <7i(t) is a non-zero constant. Therefore, (<7i(t), ■ ■ ■ ,a n (t)) algebraically integrates 
system ([5]). The initial conditions in the equation for ai(t) are as follows: oi(to) = CTi(£i(to)i ■ • ■ )£n(to)) — 
£i(to) + • ■ • + £n(to) an d o~[ (to) = (£ fc a"i)(to)- The condition &(to) ^ £j(*o) f° r an y i J guarantees 
A„ ^ 0. ' □ 

Thus, we have reduced the problem of integrability of a generic symmetric quadratic dynamical system to 
the question of solving an ordinary differential equation in h = a%. Set £ = [h, . . . , /j( n_1 )) T . This equation 
takes the form 

(15) A "C=0 

||w||=-4(n+l) 
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for the multi-index u> = (ii,...,i n ) with ||w|| = ikdegh^ k ^ and the grading degh = —4, degi — 4. For 
relations of such equations to heat equation solutions see [1] . 

Theorem 31. Each almost generic symmetric system ^ with the initial data £(io) — (^i(^o)j • • • > Cn(^o)) 
such that £,i(to) 7^ £j(^o) for any i ^ j is algebraically integrable, namely there is an ordinary differential 
equation of order n — 1, differential polynomials /12, . . . , h n ^\ in hi, and a first-order ordinary differential 
equation on h n with coefficients being differential polynomials in h\ such that for the each solution hi of 
the ordinary differential equation and each corresponding solution h n the set of functions hi, h%, ■ ■ ■ , h n with 
appropriate initial conditions algebraically integrates this system. 

Proof. According to Theorems [191 H~5l and [30l it is sufficient to prove the Theorem for n = 2 and 3 in the case 
of almost generic non-generic systems. Following the proof of theorem 1301 consider system (|6]) with = Ok 
being the basis of elementary symmetric functions. Under the conditions of the Theorem we have Ck 7^ 
for k = 1, . . . , n — 2 and c„_i = 0, thus ci = for n = 2, and ci ^ 0, c-i — for n — 3. Hence in the case 
n — 3 the function a2(t) can be expressed as a polynomial in <7i(t) and its derivatives. Thus the (n — l)-st 
equation gives a single homogeneous differential equation for o~i(t) in which the coefficient multiplying the 
highest derivative of o~i(t) is a non-zero constant. The n-th equation takes the form 

a n = coiCTn + g{ai,a x ,...,al ') 

for some constant c and some polynomial g. If c = we find a n as an integral of g{o~i, a[, . . . , cr 1 ™' ) ), else set 
o~ n = 4>i4>2, where (<f>i, fa) is a solution of the system 

4>'i = cci^i, 

±1 1 / ' ( n )\ 

</>2 = —g{°~i,vi,---,o-[ '). 

Therefore, we have described the set . . . ,a n (t)) that algebraically integrates system ©. The initial 

conditions in the equations are the same as in theorem 1301 □ 

2.2. The action of the normalizer subgroup in GL(n, C) of generic symmetric quadratic dynam- 
ical systems on the differential equation related to the problem of integrability. 

The action of the normalizer subgroup in GL(n, C) of generic symmetric quadratic dynamical systems on 
equation (|15[) is induced from the representation of GL(n,C) on dynamical systems in the following way: 
the matrix B G GL(n, C) brings £ to 77 = B£, thus the equation on h = cti(£) to the corresponding equation 

The action of S„ on equation (IT5|) is trivial. Thus the action of the normalizer subgroup in GL(n, C) of 
generic symmetric quadratic dynamical systems coincides with the action of B. For the matrix S(A, q) we 
have h = (A + nq)h. Thus equation (|15p becomes 

+ ^ = 

|M|=-4(n+l) 

for ( = (h, . . . , M Tl-1 )) T , A w = (X + nq) 1 ^^ X u , u) — (ii,...,i n ) and \uj\ = ^2 k ik- The other symmetric 
functions in this case are expressed as functions of o~k(£.) as 

= E ( n ~ k+m )* k - m r<rk- m (o<Ti(tr, 

where <ro (C) = 1 • This formulas are proved by direct calculations using the definition of . 

3. Applications. 



In this section we show how our method of algebraic integration works for finding solutions of known 
problems of quadratic dynamical systems theory. 
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3.1. One-dimensional quadratic dynamical systems. 

Each one-dimensional quadratic dynamical system has the form 

ei = < 

for some a. It is symmetric and generic for any a according to Theorem 1151 We have <7i(£) = pi(£) = £1, 
thus system §5§ for this generators as well as equation (|T5|) coincides with the dynamical system considered. 
For a 7^ the general solution to this system is 

f, 



c — at 

for some constant c. For a = the general solution is £i = c. 

The one-dimensional quadratic dynamical system is algebraically integrable by £i(t) given above. 

3.2. Two-dimensional symmetric quadratic dynamical systems. 

The general two-dimensional symmetric quadratic dynamical system has the form 

(16) ei = (« + P + 7)f? + 03 + 27)66 + 7& 

£ 2 = (a + /3 + + 08 + 2 7 )66 + 7& 

It is generic for a^O. The orbits of the action of B(X,q) on such systems are described in Lemma 121)1 
The dynamical system ([5]) for homogeneous generators <X\ — p\ and a 2 = P2 takes the form 

pi = (/3 + 2 7 )p 2 + a P2 , 

p' 2 = -(a- 2 7 )p? + (3a + 2(3) PlP2 . 

For homogeneous generators a\ = <J\ and a 2 = o 2 it takes the form 

a\ = (a + /3 + 2 7 )ct 2 _ 2acr 2 , 

^2 = 7 cr i + (a + 2/3)CTicr 2 - 
Thus, for a^O the function <ti = p\ is the solution of equation (TT51) which takes the form 

(17) /i" - (3a + 4/3 + 4 7 )M' + (a + /3)(a + 2/3 + 4 7 )/i 3 = 
and we have 

n«i * q + /? + 2 7 1 

(18) CT2 = to - ^"to* 1 ' 

In this case equation (|T4")) for hi = £7, takes the form 

£ 2 - + cr 2 = 0. 

We have A 2 = a\ — 4cr 2 . Therefore A 2 = if o\(t) satisfies the differential equation 

2ti = (a + 2/3 + 4 7 )/i 2 . 

The general two-dimensional generic symmetric quadratic dynamical system (I16|) is algebraically integrable 
by the set (<Ti(t),<7 2 (t)) where tr 1 (£) = h,(t) is a solution to (JTTJ) such that 2cr^ (t ) ^ (a + 2/3 + 4 7 )cri(t ) 2 
and (72 (t) is defined by ([T5|) . 

Let us consider special cases of equation (fTTj) . Denote Ai = —(3a + 4/3 + 4 7 ), A 2 = (a + /3)(a + 2/3 + 4 7 ), 
so that equation p7|) takes the form 

(19) h" + \ihti + \ 2 h 3 = 0. 



In the special case Ai = the general solution to (fH?|) is = fc 2 sn yyy 2 * ^ ) ^2, *J , where sn is 
the Jacobi sine, i is the square root of —1 and k\, k 2 are constants. 

In the special case A2 = the general solution to (|19[) is h(t) = -7^=- tanh (<J ^^(t + k 2 )^j , where ki, k 2 
are constants. 

In the special case A 2 = 9A2 the general solution to (fl9l) is h(t) — 6(&ii + k 2 )/Xi(kit 2 + 2k 2 t + 2), where 
ki, k 2 are constants. 



Let us now consider the case of non-generic two-dimensional symmetric quadratic dynamical system (|16p , 
that is the case a = 0. For homogeneous generators a\ = o\ and a 2 — a 2 system ([5]) takes the form 

<j[ = {P + 2 1 )al 

<t' 2 = 20<j 1 <j 2 + i<y\. 

For 6 = (0 + 2j) ^ the general solution to this system is 

<Ti = -(St + kx) , a 2 = -(5t + k 1 ) + k 2 (St + fci) -8- 

for some constants k\, k 2 . 

For (0 + 27) = the general solution to this system is 

<ti = *i, ct 2 = fc 2 e mit + -k\ 

for some constants fci, k 2 . 

Therefore in this case we get the following explicit result: 

The general two-dimensional non-generic symmetric quadratic dynamical system (|16|) is algebraically 
integrable in U — C by the set (a\(t), 02 (i)) where <j\ and 02 with k2 = are given above. 

3.3. Three-dimensional symmetric quadratic dynamical systems. 

The general three-dimensional symmetric quadratic dynamical system has the form 

£ = (a + + 7 + <% 2 + 08 + 27)^(6 + 6) + (7 + <5)(£ 2 2 + £ 3 2 ) + 2 7 66, 

(20) & = (a + + 7 + 5)e 2 + (0 + 2 7 )6(6 + &) + (7 + ml + & + 27^6, 

& = (a + /3 + 7 + 5)£ 3 2 + 08 + 27)6(6 + 6) + (7 + <5)(C 2 + Cl) + 2766- 

It is almost generic for a + 36 ^ and generic for a^O and a + 3<5 ^ 0. 
Equation (fT4|) takes the form 

e 3 - /m£ 2 + ^ - h 3 = 

with A 3 = -27hj + 18hih 2 h 3 - Ah\h 3 - Ah\ + h\h\. 

The dynamical system ([5]) for homogeneous generators ai — <J\, a 2 — <r 2l a 3 — a 3 takes the form 

a[ = (a + + 3 7 + 36)aj - 2(a + 3<5)<7 2 , 

4 = 2(7 + S)af + (a + 20 - i6)a 1 a 2 - 3aa 3 , 

O3 = (7 + 5)cr2Cr 2 — 2(5cT2 + (« + 3/3)(T3C7i. 

The general three-dimensional generic symmetric quadratic dynamical system (|20[) is algebraically inte- 
grable by the set (a\(t), cr 2 (t), a 3 (t)) where ci is a solution of equation 

(21) A1/1'" + A 2 /i"/i + \ 3 {tif + Xih'h 2 + X 5 h 4 = 
with 

Ai = a + 3i5, 

A 2 = -Ai(4a + 7/3 + 6 7 + 2<5), 

A 3 = -(3a 2 + 6S 2 ) - 2(2a0 + 3a 7 + AaS + 6/3(5 + 9 7 <5), 

A 4 = 6a 3 + 2(lla 2 /3 + 9^ 2 a + 15a 2 7 + 13a 2 <5 + 12a<5 2 + 27/3 2 6 + 1805 2 ) + 36(a0j + 2a0S + 2ajS + 30jS), 

A 5 = -(a + 30 + 9 7 + 35)(a 3 + 3a 2 + 2a0 2 + a 2 j + 3a 2 5 + 60 2 S + 8a0S) 

and 

-2(a + 3S)a 2 = a[ - (a + + 3 7 + 3<5)<j 2 , 

6a(a + 35)a 3 = o'{ - (3a + A0 + 67 + 2S)a[ai + (a 2 + 20 2 + 3a0 + 7a 7 + 3aS + 6/3 7 + 20S)af, 

where the initial conditions satisfy 27ct3(£o) 2 7^ ^8ai(to)a 2 (to)(r 3 (to)— iai(to) 3 cr 3 (to)— Acr 2 (toY+(Ji(to) 2 (j 2 (to) 2 . 
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Let us now consider the case of almost generic non-generic three-dimensional symmetric quadratic dy- 
namical system (|20p . that is the case a = 0, 5 0. For homogeneous generators a\ = a±, ai — o~2j 03 = 03 
system §5$ takes the form 

o-; = (/? + 3 7 + 3<5)<j2 _ efca, 
(7 2 = 2(7 + S)af + 2(0 - 25)a 1 a 2 , 

03 = (7 + <5)(7 2 (72 - 2(5(7 2 + 3/3<7 3 C7l. 

The almost generic non-generic three-dimensional symmetric quadratic dynamical system (|20[) is alge- 
braically integrable by the set (a±(t), 02 (i), 03 (i)) where <ti is a solution of equation 

h"(t) - 2(2/3 + 3 7 + 6)h(t)h'(t) + 2/3(fi + 3 7 + <5)/i(t) 3 = 0, 

we have 02 — ^ +3 ^ +3 ^ u\ — gyOi, and for this a± the function (73 is a solution of equation 

I8&73 - 54/3(5cr 3 (7i + (cri) 2 - (2/3 + 3 7 + 3(5)cr 2 c4 + /3(/3 + 3 7 + 3£)oi = 0. 

The equation on initial conditions is the same as in the previous case. 

3.4. Lotka-Volterra type system. 

Let us consider the dynamical system 

(22) e* = &fe6-2&Y fc = l,2,3. 

It is symmetric and has the form (|20[) for a = —2, /3 = 1, 7 = 0, S = 0, thus it is generic. 

This system is algebraically integrable by the set of functions (01 (t), o~2(t), o~3{t)) where 01 is a solution 
of equation (see (f2Tj) ) 

(23) + fi/'/i + 2(/i') 2 - 2h'h 2 = 

and 4(7 2 = 0"i + cr 2 , -64cr 3 = cr^ + 2a[ai. 

S. Kovalevskaya introduced system (j2"2")l in a letter to G. Mittag-LefHer. She obtained the result that this 
system is integrable in elliptic functions. It has the solution £1 = (fnsn(t))', £2 = (lncn(i))', £3 = (In dn(t))' . 
We have o\ = £1 + £2 + £3, therefore the elliptic function (ln((sn(i) 2 )'))' is a solution of ([23]) . 

For dimensions n)3a generalization of system (|22[) is a Lotka-Volterra type system (see [3]) 

(24) 4=&^^- 2 ^)' k = l,...,n. 

This system is symmetric and generic for any n, thus our general method gives the answer to the question 
of algebraical integrability for this system. 

Example 32. For n = 4 equation (fT5|) takes the form 

ft"" _ ft"'/! + 5h"ti - 4h"/i 2 - 8(h') 2 h + ih'h 3 = 

and the differential polynomials are 



(2-,) <r2 = jK+<7?), ,73^1(^ + 2(7^), ^^K'+<t' 1 Vi + 2K) 2 -2«tW) 



Lemma 33. The system (|24|) /or n = 2 /ias onZy one quadratic integral x\Xi- 
For n — 3 it has two independent quadratic integrals 

ax\Xi + bx\X 3 + CX2X3, a + b + c = 0. 

For n = 4 it has two independent quadratic integrals 

a(x\X2 + X3X4) + b(xix 3 + X2X4) + c(x 2 x 3 + X1X4), a + b + c = 0. 

For n ^ 5 the system (|24p has no quadratic integrals. 



Proof. System (J2U) is such that CP{x\, . . . , x n -i, 0) = (£P(x\, . . . , ar n _i, x n )) \ x =Q . Therefore for an inte- 
gral P(xi , . . • , i„_i, x n ) of system (|24|) for some n the polynomial P(xi , . . . , x n —i, 0) is an integral of system 
([24)1 for n - 1. 

One can check directly that the lemma holds for n = 4 and n = 5, all other cases becoming corollaries of 
the last statement. □ 

Another approach to integration of system ([24]) (see 4 ) is a relation of this system to a generalization of 
the Eulcr top 

f?/b = n>'' k = l,...,n, 
using the transform £fc = Yli^k Vi- One gets n — 1 functionally independant rational integrals of the form 



I C C \^ n ~2 n 



it 



fc- 

fc=l 



Remark 34. TTie system 

(26) £[• - 6 ^ & " m ^ > * = 1, 2, 3, 

/or m ^ 2 /ias no quadratic integrals. 

3.5. The Darboux-Halphen system. 

Let us consider the classical Darboux-Halphen system 

& = 66 - 66 _ 66 j 
£2 = 66 - 66 - 66) 
6 = 66 - 66 - 66- 

It is symmetric and has the form ([20]) for a = 2, (3 = —2, 7 = ^,0" = — |, thus it is generic. 

The Darboux-Halphen system is algebraically integrable by the set of functions (o"i(i), o"2(i), o"3(i)) where 
o-\ is a solution of equation (see (|2ip) 

ft'" + 4/ l "fr-6(ft') 2 = 

and o"2 = — cr'x, 60-3 = o'{. Therefore we have o~\ = — hy, where y(t) is a solution of the famous Chazy-3 
equation 

y"> -2y"y + 3(y>) 2 =0. 
Let us consider a generalization of the Darboux-Halphen system, namely 

£ = a(6&-6£i + 
£ 2 = a(66- 66-66) + ^, 
6 = a(66-66-66) + 66 2 - 

It is symmetric and has the form (|20p for a = 2a + 0, /3 = —2a, 7 = |a, 5 = —\a. thus it is generic for 
2a ^ —6, a ^ — 2o. It coincides with the classical Darboux-Halphen system for a = 1, b = 0. 

According to (f2~Tj) . for a 7^ for this system we have 2(a — b)a\ — —y, where y(t) is a solution of 

(27) y / "-2y y " + 3(y>) 2 +cW-y 2 ) 2 =0, 



i(a + 2b)(a-b) 



This non-linear ordinary differential equation coincides with the Chazy-3 equation for c = and with the 
Chazy-12 equation for c = fc ^ 36 , fc € N, fc > 1, fc 7^ 6. These equations are examples of known third-order 
differential equations with the Painleve property. 
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3.6. An almost generic classical system. 

Let us consider the classical three-dimensional dynamical system in Lax form 

(28) L' = [[K,L k ],L], 



where L = — £1 £3 and K is a 3 x 3 constant matrix. This system is quadratic for k = 1. In 




d mi 777,2 

this case for the equality of the right and the left part in (|28]l it is necessary that K = \ —mi d 7773 

-7772 — ITI3 d 

for some d, mi, 7772, 7773. Moreover, the resulting system does not depend on d, thus we can put d = 0. The 
system is symmetric if and only if 777,1 = m2 = "^3 = m. Therefore any symmetric quadratic dynamical 
system of the form (1281) up to a rescaling £j — > *l coincides with the symmetric quadratic dynamical system 

£i=£i(6 + £3)-fe 2 +£ 3 2 ), 

(29) £ 2 = 6(£i + £3)-(£ 2 + £3 2 ), 

£3 = 6(£i + 6)-(£ 2 + £2 2 )- 

This system has the form (|20|) for a = 0, /3 = 1, 7 = 0, 5 = — 1. It is not generic, but almost generic. 

From the form (j2"5|) we get integrals of the system of the form tr(L k ). We have tr(L 2k+1 ) = and 
tr(L 2k ) — Cfc(£i + £! + £i) fe f° r constant with ci = —2, thus this method gives only one integral equal 
P2(0- Let us use our method to algebracally integrate this system. 

For system ([29)) the funtion a\ is a solution of equation 

(30) h" - 2h'h = 

and (T2 = + 2<Ti). The expression a\ — 2<T2 = P2 is an integral of the system. For a solution o\ of (|3U)) 
the function 03 is a solution of 

(31) a' 3 = 3a 1 a 3 + ^(a[-a 2 )(a[+2a 2 ). 

The general solution of (|30[) is tri = ctan(c7; + &). We get 02 = — ^- ^2 — cos (ct+b) 2 ) • ^ ne corre sponding 
general solution to (|31~|) is 

c 3 sin(ct + 6)(2 cos(ct + b) 2 - 5) + k 
° Z ~ 54cos(c< + b) 3 

for a constant k. 

System (|2"9"|) is algebraically integrable by the set (a\(t) 1 a2it) , a^t)) where cr\, d2 and a$ are given above 
and 2c 6 + k 2 ^ 0. 
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